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Abstract —With huge unlicensed bandwidth available in most 
parts of the world, millimeter wave (mmWave) communications 
in the 60 GHz band has been considered as one of the most 
promising candidates to support multi-gigabit wireless services. 
Due to high propagation loss of mm Wave channels, beamforming 
is likely to become adopted as an essential technique. Conse¬ 
quently, transmission in 60 GHz band is inherently directional. 
Directivity enables concurrent transmissions (spatial reuse), 
which can be fully exploited to improve network capacity. In this 
paper, we propose a multiple paths multi-hop scheduling scheme, 
termed MPMH, for mm Wave wireless personal area networks, 
where the traffic across links of low channel quality is transmitted 
throngh mnitiple paths of multiple hops to unleash the potential 
of spatial reuse. We formulate the problem of multiple paths 
multi-hop scheduling as a mixed Integer linear program (MILP), 
which is generally NP-hard. To enable the implementation of the 
multiple paths multi-hop transmission in practice, we propose 
a heuristic scheme including path selection, traffic distribution, 
and multiple paths multi-hop scheduling to efficiently solve the 
formulated problem. Finally, through extensive simulations, we 
demonstrate MPMH achieves near-optimal network performance 
in terms of transmission delay and throughpnt, and enhances 
the network performance significantly compared with existing 
protocols. 


I. Introduction 

Recently, millimeter wave (mmWave) communications in 
the 60 GHz band has attracted considerable interest from 
academia, industry, and standards bodies. Due to its huge 
unlicensed bandwidth (i.e., up to 7 GHz in the USA), 60 
GHz communications is able to support multi-gigabit wireless 
services, such as high-speed data transfer between devices 
(e.g., cameras, pads, and personal computers), wireless gigabit 
ethernet, wireless gaming, and real-time streaming of both 
the compressed and uncompressed high definition television. 
Meanwhile, rapid progress in 60-GHz mm-wave circuits. 
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including on-chip and in-package antennas, radio frequency 
(RF) power amplifiers (PAs), low-noise amplifiers (LNAs), 
voltage-controlled oscillators (VCOs), mixers, and analog-to- 
digital converters (ADCs) has accelerated popularization of 
wireless products and services in the 60 GHz band |[T|-|I3- In 
terms of standardization, several standards have been defined 
for indoor wireless personal area networks (WPAN) and 
wireless local area networks (WLAN), for example, ECMA- 
387 El , IEEE 802.15.3c E), and IEEE 802.1 lad Ig). Due 
to high carrier frequency, mmWave communication systems 
are fundamentally different from other existing communication 
systems using lower carrier frequencies (e.g., from 900 MHz 
to 5 GHz). Therefore, network architectures and protocols, 
especially medium access control (MAC) mechanism, require 
new thinking to fully unleash the potential of mmWave com¬ 
munications. 

MmWave communications in the 60 GHz band suffers from 
high propagation loss. The free space propagation loss scales 
as the square of the wavelength, which indicates that the free 
space propagation loss at 60 GHz band is 28 decibels (dB) 
more than that at 2.4 GHz Q. Thus, 60 GHz communica¬ 
tions is mainly used for short-range indoor communications. 
To combat severe channel attenuation, high gain directional 
antennas are utilized at both the transmitter and receiver, where 
beamforming has been adopted as an essential technique to 
search for the best antenna weight vectors (AWV) ISl- lfTOll . 

With directional transmission, the third party nodes cannot 
perform carrier sense as in IEEE 802.11 to avoid contention 
with current transmission, which is know as the deafness 
problem im. On the other hand, there is less interference 
between links, and thus concurrent transmission (spatial reuse) 
can be exploited to greatly increase the network capacity. 
However, on one hand, for flows across links of low channel 
quality or with significantly higher traffic demand, more 
network resources (e.g., time slots) are needed. In this case, 
these extra resources cannot be utilized efficiently due to less 
spatial reuse. Consequently, system performance is degraded 
significantly. Thus, more efficient transmission mechanisms 
are needed to address the performance degradation when 
there are flows across links of low channel quality or with 
significantly higher traffic demand. Eor example, in na, a 
traffic flow of a long hop is transmitted over multiple short 
hops to improve flow throughput. Similarly, in D-CoopMAC 
Ha, the direct long link is replaced by a two-hop link if a 
relay that has higher rate links with both source and destination 
exists. On the other hand, for flows supporting applications of 
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high throughput, such as high-definition television (HDTV), 
much more time slots are needed to satisfy their throughput re¬ 
quirements in TDMA-based protocols such as IEEE 802.15.3c. 
Consequently, there are not enough time slots allocated for 
other flows, which may lead to serious unfairness in WPANs. 
Besides, the number of flows scheduled each time in these 
protocols will be very small. Therefore, fully exploiting spatial 
reuse to improve the throughput of these flows while achieving 
high system performance is an important and challenging 
issue. 

Recently, the hybrid beamforming structure has been pro¬ 
posed to obtain the multiplexing gain of Multiple-Input 
Multiple-Output (MIMO) and also provide high beamforming 
gain to overcome high propagation loss in mmWave bands 
M- In the hybrid beamforming structure, multiple RE links 
can be established via the design of the digital precoder and 
the analog beamformer to reap the spatial multiplexing gain 
of MIMO HSl, Hb). To improve the throughput of flows of 
low channel quality or with high traffic demand, transmitting 
these flows through multiple paths concurrently similar to the 
multiple RE links in the hybrid beamforming structure will be 
an effective way, and should be considered in the design of 
scheduling schemes for mmWave communications. 

In this paper, motivated by the spatial multiplexing structure 
in hybrid beamforming, we propose a novel multiple paths 
multi-hop scheduling scheme (MPMH) to boost the spatial 
reuse in mmWave WPANs. By transmitting the traffic of flows 
with low channel quality on the direct paths or high traffic 
demand through multiple multi-hop paths, the potential of spa¬ 
tial reuse is unleashed, and concurrent transmissions on these 
paths are fully exploited to enhance the system performance 
in terms of flow and network throughput. The contributions of 
this paper are three-fold, which are summarized as follows. 

• We formulate the multiple paths multi-hop scheduling 
problem into a mixed integer linear program (MILP), 
i.e., to minimize the number of time slots by multiple 
paths multi-hop transmissions with the traffic demand of 
all flows accommodated. Concurrent transmissions, i.e., 
spatial reuse, are explicitly considered under the signal to 
interference plus noise ratio (SINR) interference model in 
this formulated problem. 

• We propose an efficient and practical scheme to solve 
the formulated NP-hard problem by three heuristic algo¬ 
rithms of path selection, traffic distribution, and trans¬ 
mission scheduling. In the transmission scheduling algo¬ 
rithm, concurrent transmissions are enabled if the SINR 
of each link is able to support its transmission rate. 

• Extensive simulations are carried out to demonstrate 
the near-optimal network performance of MPMH, and 
superior performance in terms of delay and throughput 
compared with other existing protocols. Performance 
under different maximum number of hops on paths is 
also investigated to provide guidelines for the choice of 
this parameter in practice. 

The rest of this paper is organized as follows. We dis¬ 
cuss related work on concurrent transmission scheduling for 
mmWave WPANs in section |II] In section HID we present the 


system model and illustrate our basic idea by an example. 
We formulate the optimal multiple paths multi-hop scheduling 
problem as an MILP in section |IV] The MPMH scheme 
is illustrated in section |V] In section IVIl we evaluate the 
performance of MPMH under various traffic patterns and 
simulation parameters. Einally, we conclude this paper in 
section Ivnl 

H. Related Work 

Transmission scheduling for mmWave communications in 
the 60 GHz band has been investigated in the literature im- 
Cl, El-Gl, ES-EtI. Since ECMA-387 H and IEEE 
802.15.3c 111 adopted time division multiple access (TDMA), 
some work is also based on TDMA ifTTll . ifTSl . E2l . EHj Efill . 
EH. In two protocols based on IEEE 802.15.3c, multiple links 
are scheduled to communicate in the same slot if the multi-user 
interference (MUI) is below a specific threshold JTtI . ESI. Cai 
et al. E 2 I introduced the concept of exclusive region (ER) to 
enable concurrent transmissions, and derived the ER condi¬ 
tions that concurrent transmissions always outperform TDMA. 
Qiao et al. E3l proposed a concurrent transmission scheduling 
algorithm for an indoor IEEE 802.15.3c WPAN, where non¬ 
interfering and interfering links are scheduled to transmit 
concurrently to maximize the number of flows with the quality 
of service requirement of each flow satisfied. Eurthermore, 
multi-hop concurrent transmissions (MHCT) are enabled to 
address the link outage problem (blockage) and to combat 
huge path loss to improve flow throughput El. MHCT, 
however, only transmits traffic of flows through multiple hops 
on one path, and does not exploit the concurrent transmissions 
on multiple paths to improve flow throughput and network 
throughput. In IEEE 802.15.3c, during the random access 
period, the piconet controller is in the omni-directional mode 
to solve the deafness problem, which may not be feasible for 
mmWave systems of the multi-gigabit domain with highly 
directional transmissions. Eor bursty traffic patterns such as 
the Interrupted Poisson Process (IPP), TDMA based protocols 
may result in over-allocated medium access time for some 
users while resulting in under-allocated medium access time 
for others. 

There is some other work based on a central controller 
to coordinate the transmissions in WPANs ED, E3, El, 
El, ED. Gong et al. ED proposed a directional carrier 
sense multiple access with collision avoidance (CSMA/CA) 
protocol, which mainly focuses on solving the deafness prob¬ 
lem. It exploits virtual carrier sensing, and depends on the 
piconet coordinator (PNC) to distribute the network allocation 
vector (NAV) information. However, it does not exploit the 
spatial reuse fully and also does not consider multiple paths 
multi-hop transmissions. In the multi-hop relay directional 
MAC protocol (MRDMAC), if a wireless terminal (WT) is 
lost, the access point (AP) will choose a WT among the 
live WTs to act as a relay to the lost node EOI . However, 
MRDMAC does not fully exploit spatial reuse since most 
transmissions go through the piconet coordinator (PNC). Chen 
et al. E 4 I proposed a spatial reuse strategy to schedule two 
different service periods (SPs) to overlap with each other 
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for an IEEE 802.11 ad WPAN. Recently, Son et al. HI] 
proposed a frame based directional MAC protocol (EDMAC). 
The high efficiency of EDMAC is achieved by amortizing the 
scheduling overhead over multiple concurrent transmissions 
in a row. The core of EDMAC is the Greedy Coloring algo¬ 
rithm, which fully exploits spatial reuse and greatly improves 
the network throughput compared with MRDMAC ESI and 
memory-guided directional MAC (MDMAC) 11211 . EDMAC 
also has a good fairness performance and low complexity. 
EDMAC, however, does not consider multiple paths multi-hop 
transmission to improve flow throughput of links with poor 
quality. Chen et al. ini proposed a directional cooperative 
MAC protocol, termed D-CoopMAC, to coordinate the uplink 
channel access among stations in an IEEE 802.1 lad WLAN. 
In D-CoopMAC, the multirate capability of links is exploited 
to select a relay station for the direct link; when the two-hop 
link outperforms the direct link, the latter will be replaced 
by the former. In D-CoopMAC, however, spatial reuse is not 
considered since most transmissions go through the access 
point (AP). Most of the above work ignores the negative effect 
of links with poor channel quality on network performance, 
such as delay, throughput, and the number of flows scheduled 
each time. The work above does not boost the potential of 
spatial reuse via multiple paths multi-hop transmissions. 

There is also some related work on maximum independent 
set and protocol model based scheduling. Eor protocol model 
based scheduling, interference is modeled as an interference 
graph, where each vertex represents a link in the wireless 
network EII-EQI- Two vertices form an edge if these two 
links cannot be scheduled for concurrent transmissions. This 
simple interference model does not take the unique features of 
mm Wave links, e.g., directivity, into consideration. Xu et al. 
ED proposed a constant approximation algorithm for one- 
slot link scheduling in arbitrary networks under the SINR 
interference model, where a transmission is successful if the 
received SINR is more than some threshold. However, it 
does not consider the directivity of mmWave links, and the 
interference is only related to the distance between nodes in a 
two-dimensional Euclidean space. This lack of consideration is 
unreasonable for mm Wave WPANs, especially since the SINR 
thresholds for all links were kept the same. To the best of our 
knowledge, we are the first to exploit multiple paths multi¬ 
hop transmissions for boosting the spatial reuse gain under 
the SINR interference model. 

III. System Overview 

A. System Model 

We consider a typical indoor WPAN system composed 
of several wireless nodes (WNs) and a single piconet con¬ 
troller (PNC) El. The PNC synchronizes the clocks of other 
nodes and schedules the medium access of all the nodes 
to accommodate their traffic demands. The WNs and PNC 
have electronically steerable directional antennas to support 
directional transmissions between WNs or between the WN 
and the PNC by beamforming. The system is partitioned 
into non-overlapping time slots of equal length, and runs 
a bootstrapping program EDj by which each device knows 


the up-to-date network topology and the location information 
of other devices. With this information, the beamforming 
between nodes can be completed in a short time, and each 
node can direct its antenna towards other nodes. 

We assume there are V flows with traffic transmission 
demand in the network. Eor flow v, we denote its traffic 
demand as and numerically dy is equal to the number 
of packets to be transmitted. Eor 60 GHz wireless channels, 
non-line-of-sight (NLOS) transmissions suffer from higher 
attenuation than line-of-sight (LOS) transmissions m, El- 
In Ref. ED, the path loss exponent in the LOS hall is 2.17, 
while the path loss exponent in the NLOS hall is 3.01. Thus, 
if the distance between the transmitter and the receiver is 10 
m, the gap between the path loss of LOS hall and NLOS hall 
is about 10 dB. In a power-limited regime, a 10 dB power loss 
requires a 10-fold reduction of transmission rate to maintain 
the same reliability. On the other hand, NLOS transmission in 
the 60 GHz band also suffers from a shortage of multipath, and 
restricting to the LOS path can maximize the power efficiency 
since the LOS path is strongest ll^ . To achieve high data rate 
transmission and maximize the power efficiency, we consider 
the directional LOS transmission case in this paper. Eor each 
directional link i, we denote its sender and receiver as Si and 
Vi respectively. Then according to the path loss model ES^ the 
received signal power, denoted by Py (mW), can be calculated 
as 

Py = k0Ptl7j,, ( 1 ) 

where Pt (mW) is the transmission power, ko = 
is the constant scaling factor corresponding to the reference 
path loss PL{do) (dB) with do equal to 1 m, (m) is the 
distance between node Si and node r^, and 7 is the path loss 
exponent ll23l . In this paper, we assume fixed transmission 
power. 

We denote the transmission rate of the direct link of flow 
V as Cy, and numerically c„ is equal to the number of packets 
the link can transmit in one time slot. The transmission 
rates of links are obtained by a channel transmission rate 
measurement procedure ES. In this procedure, the sender 
of each flow transmits measurement packets to the receiver 
first. After measuring the signal to noise ratio (SNR) of these 
packets, the receiver obtains the achievable transmission rate, 
and appropriate modulation and coding scheme (MCS) by 
the correspondence table about the SNR and MCS. Under 
low user mobility, the procedure is usually executed, and the 
transmission rates of links are updated periodically. 

Directional transmissions enable less interference between 
links, and concurrent transmissions of links can be supported 
to improve network capacity. Due to the limited range for 
mmWave WPANs, the interference between links cannot be 
neglected E2l . Thus, we adopt the SINR model for concurrent 
transmission El, EH- Eor link i and link j, the received 
power from Si to Vj can be calculated as 


where fsi,yj indicates whether Si and rj direct their beams 
towards each other. If Si and Vj direct their beams towards each 
other, fsi.Tj = 1; otherwise, fsi,yj = 0. Then, the received 
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SINR at rj, denoted by SINRj, can be calculated as 


SINRj 


WNo+pY.fsur,koPtl7j/ 

i¥^3 


(3) 


where p is the multi-user interference (MUI) factor related to 
the cross correlation of signals from different links, W (Hz) 
is the bandwidth, and Nq (mW/Hz) is the one-sided power 
spectra density of white Gaussian noise m- For each link 
j, we denote the minimum SINR to support its transmission 
rate Cj as MS{cj). Therefore, concurrent transmissions are 
supported if the SINR of each link j is larger than or equal 
to MS{cj). 

In MPMH, we adopt a frame based scheduling scheme, 
which is shown in Fig. [Ua). Node A is the PNC, and other 
nodes are WNs. In MPMH, time is divided into a sequence 
of non-overlapping frames ED. Each frame consists of a 
scheduling phase and a transmission phase. In the scheduling 
phase, after all of the WNs steer their antennas towards the 
PNC, the PNC polls the WNs successively for their traffic 
demand in time tpoii- Then the PNC computes a schedule to 
accommodate the traffic demand in time tsch- Finally, the PNC 
pushes the schedule to WNs in time tpush- In the transmission 
phase, all devices communicate with each other following the 
schedule until their traffic demand is cleared. Spatial reuse 
gain can be boosted via multiple paths multi-hop transmission 
in the transmission phase to improve flow as well as network 
throughput, which depends on the scheduling solution. 


B. Problem Overview 

For flows with low channel quality on their direct paths or 
high traffic demand, to improve throughput, their traffic should 
be transmitted through multiple paths to unleash the potential 
of spatial reuse of hops on these paths. Meanwhile, the traffic 
demand of flows should be accommodated with a minimum 
number of time slots in the transmission phase to maximize 
the transmission efficiency. 

Now, we present an example to illustrate the operation of 
MPMH and our basic idea. We assume a WPAN of six devices 
in the network, denoted as A, B, C, D, E, and F, as shown 
in Fig.lHb). The numbers above the directional edges between 
nodes represent their transmission rates, and numerically are 
equal to the numbers of packets these links can transmit in one 
time slot. If there are 18 packets to be transmitted from A to 
B, then with a transmission rate of I packet per time slot, 18 
slots are needed in the transmission phase to clear this traffic 
demand. By MPMH, we select three paths from A to B, i.e., 
the direct path from A ^ B (path 1), the path of A —>• C —>• 
E ^ B (path 2), and the path of A^D^F^B (path 
3). Then we distribute 3 packets to path 1, 9 packets to path 
2, and 6 packets to path 3. Afterwards, MPMH computes a 
schedule as illustrated in Fig. [11^), where each colored box 
represents a time slot. The schedule has 5 pairings, and in the 
third pairing of the schedule, links A —3 B, C ^ E, and D 
F transmit concurrently for 3 time slots. This schedule clears 
the packets from A to B in 9 time slots, and by transmitting 
this flow on three paths, concurrent transmissions of the three 


links in the third pairing are enabled to improve the efficiency 
of transmission significantly. 

IV. Problem Formulation 

In this section, we formulate the problem of optimal mul¬ 
tiple paths multi-hop scheduling into a mixed integer linear 
program (MILP) based on the problem formulation in FDMAC 

ini. 


A. Problem Formulation and Analysis 

From traffic demand polling, we assume there are V flows 
to be scheduled by the PNC. For flows across links with a 
low transmission rate or with a high traffic demand, more 
time slots are needed to accommodate their traffic demand. 
Consequently, these extra time slots cannot be used for spatial 
reuse, which degrades system performance significantly. Thus, 
flows with lower transmission rate on the direct paths or with 
higher traffic demand have higher priority to be transmitted 
through multiple paths for better utilization of spatial reuse. 
In the following, we propose a criterion to decide whether a 
flow needs to be transmitted through multiple paths. 

For each flow v, we define its traffic demand intensity as 
arrived traffic demand averaged over a relatively long time, 
which is denoted by By. The transmission rate of its direct 
path is denoted by Cy. Then if flow v satisfies 


Cy / Dy 

avg(c„/D„) 

uev 


(4) 


then traffic of flow v will be transmitted through multiple 
paths, e is defined to control the number of flows transmitted 
through multiple paths. With a larger e, there may be more 
flows transmitted through multiple paths; otherwise, fewer 
flows will be transmitted through multiple paths. Specially, 
if Cy is equal to 0, i.e., the direct path of flow v is blocked, 
flow V will be transmitted through multiple paths. 

For the fairness among flows, if there is no multi-path multi¬ 
hop transmission, flows across direct paths of low channel 
quality or with high traffic demand will occupy a large number 
of time slots exclusively. With multi-path multi-hop transmis¬ 
sion enabled, more flows are able to share the time resources 
more fully by the concurrent transmissions (spatial reuse). 
Thus, the fairness performance without multi-path multi-hop 
transmission will be worse than that with multi-path multi¬ 
hop transmission. Therefore, our scheme improves the fairness 
among flows by the multi-path multi-hop transmission. 

We denote the number of paths of the flow as My. We 
also denote the number of hops of the path of flow v as 
Hyp. For flow V, the traffic demand distributed to the path 
is denoted as dyp. We denote the hop link of the path 
of flow V as (u,p, i), and its transmission rate by Cypi. For link 
{v,p,i) and {u,q,j), we define an indicator variable lypi^uqj- 
Ivpi.uqj is equal to 1 if and {u,q,j) are adjacent, i.e., 

have common vertices; otherwise, lypi^uqj is equal to 0. 

If there are K pairings in the schedule to accommodate the 
traffic demand of flows, we denote the number of time slots 
of the pairing as 5^ ifTTl . We also define Qy ^ to indicate 
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(a) A MPMH frame 

Fig. 1. An example of MPMH operation in a WPAN of six nodes. 


(b) Transmission path illustration of MPMH 


whether link {v,p,i) is scheduled in the pairing. If link 
{v,p,i) is scheduled in the pairing, is equal to 1; 
otherwise, is equal to 0. To optimize system performance, 
the traffic demand of flows should be cleared in the shortest 
time HD. Thus, the problem of optimal multiple paths multi¬ 
hop scheduling (PI) is formulated as follows. 


K 


k=l 


(5) 


s. t. 


^ = 1, if dvr, > 0 k i < Hy 




k ^ ) {0, 1}, if dyp 0 & i ^ Hyp^ 


^vpi S t |Q|^ otherwise; 


yv,p,i,k (7) 


K > 

k=l = 0 


M„ 

dyp = dy ] V V 


Cvpi 


, if dyp ^ 0 & z ^ Hyp 1 w 

^ ^ V z;,p, i 

otherwise; 

( 8 ) 

(9) 


E o.vp^ < 1; V v,p,k 


( 10 ) 


■ -I- 1 \f T - ■ — 1 * 

'-^vpi ' ^uqj — ^vpi,uqj — 


for all A:, any two links (v,p, z), (zz, q^j) 


( 11 ) 


X! ^vpi — S ^ l^i 


^vpi — ^vp{i-\-l) 

k—1 k—1 


\fv,p,i = l^ {Hvp-l), K* = l^K 


( 12 ) 


koPtl7~' . r -C 

^vpi > 'vpt 


V Mu ^uq 

WNo+pY: E E 

14 = 1 q=l j = l 

> MS{cvpi) X M v,p,i,k 


(13) 


These constraints are explained as follows. 

• Constraint (|6]l indicates for each link {v,p,i), if there is 
traffic distributed to it, then it should be scheduled once 
in one pairing of the frame. 

• Constraint (|7]i indicates for each link {v,p,i), if there 
is traffic distributed to it, then is a binary variable; 
otherwise, ^ is equal to 0. 

• Constraint (l8]l indicates for each link {y,p,i), if there 
is traffic distributed to it, then this traffic should be 
accommodated in the frame. 

• Constraint indicates for flow v, the sum of traffic 
distributed in all paths should be equal to its traffic dy 

• Constraint (fTOl) indicates links in the same path cannot be 
scheduled in the same pairing due to the inherent order 
of transmission in each path. Preceding hops should be 
scheduled ahead of hops behind on the path since nodes 
behind on the path are able to relay the packets after 
receiving the packets from preceding nodes. 

« Constraint (fTTI) indicates due to the half-duplex assump¬ 
tion, adjacent links cannot be scheduled in the same 
pairing. 

• Constraint (O indicates due to the inherent order of 
transmission in each path, the hop link of the 
path of flow V should be scheduled ahead of the (i -f 1)*^ 
hop link. Constraint (fT2l) represents a group of constraints 
since K* varies from 1 to K. 

• Constraint (O indicates to enable concurrent transmis¬ 
sions, the SINR of each link in the same pairing should 
be able to support its transmission rate. 


B. Problem Reformulation 

We can observe that in the constraints (|6l), (|7]), and ([8]l of 
problem (PI), the condition “if dyp > 0” has the variable 
dyp. This problem form is intractable. If we select the paths 
for each flow that needs multiple paths multi-hop transmission 
by a heuristic path selection algorithm, these paths will have 
traffic, and their conditions “if dyp > 0” will be met. In this 
case, we can remove “if dyp > 0” from constraints dSI, ©, 
and ®, and problem (PI) becomes a mixed integer nonlinear 
program (MINLP). However, it is still difficult to obtain the 
optimal solution since constraints ([8]) and (fTSl) are nonlinear. 
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If the nonlinear terms of problem (PI) can be linearized, 
problem (PI) will become a standard mixed integer linear 
program (MILP), which can be solved by some existing so¬ 
phisticated algorithms, such as the branch-and-bound method 
061 . For the second order terms in constraints ® and (fT3]) . we 
adopt a relaxation technique, the Reformulation-Linearization 
Technique (RLT) HU, [361 to linearize constraints ([8]l and 
(El. The RLT technique produces tight linear programming 
relaxations for an underlying nonlinear and non-convex poly¬ 
nomial programming problem. In the procedure, a variable 
substitution is applied for each nonlinear term to linearize the 
objective function and the constraints of this problem. Besides, 
nonlinear implied constraints for each substitution variable are 
generated by taking the products of bounding terms of the 
decision variables, up to a suitable order. Specially, we present 
the RLT procedures for (H and (fTl in the following. 

For constraint (H, we define a substitution variable sj) ^ = 




vpi 

Since the traffic of one link should be accom- 


6^ 

modated in one pairing, the number of time slots of each 
pairing, S^, is bounded as 0 < 6^ < d, where d = 
max{ |for all v,p,i]. We also know 0 < < 1. 

Then we can obtain the RLT bound-factor product constraints 
for as 


{[b'" - 0] • - 0 ]}ls > 0 

{[d- s’"] ■ [a^ i - 0]}ls > 0 


{[b'^ - 0] • [1 - a„pj]}Ls > 0 


, {[d-S ] ■ [1 - <pJ}ls > 0 


for all v,p,i,k. (14) 


{■}ls represents a linearization step under 
By substituting sj ^ = S’" ■ aj ,, we obtain 


= S'" 


s > 0 

^vpt — 

i’" - s’" 

vpi vpi — 

f - ^ 0 

d-s’"-d-a’;^. 


for all v,p,i,k. (15) 


> 0 


vpi — 

For constraint ( fT3l l, we first convert it to 


{koPtl 


-7 

vpi j ' vpi 


-MS{c,p,)WNo)xatp, > 


V Mu Hug 

MS{ 

^vpi )pEEE^ 

u—1q—1j—1 


iqj ,1 vpi Vpi Uqj 


^uai a ■ r 

uqj ^uqj »' vp% 


(16) 


For the second order term we define ut’^pi^^qj = 

^vpi^uqj '^6e substitution variable. Since 0 < a’^p^ < 1 and 
0 < < 1, the RLT bound-factor product constraints for 

^vpi,uqj are 


^vpi,uqj — ^ 

> f) 

^vpi ^ vpi,uqj — 

r,^ _ > f) 

^uqj vpi,uqj — ^ 

1 - aL,- - oL, + 


(17) 


vpi 


^uqj 


vpi,uqj — 


> 0 


After substituting s’^p^ and uj’^pi^^gj into constraint ([8]l and 
constraint (El, we obtain a mixed integer linear program 
(MILP) relaxation (P2) as 


K 

min^(5'= 

k^l 

S. t. 


K 


k=l 



> 



if i < iT. 




0, otherwise; 


V v,p,i 


(18) 


(19) 


{koPtlsjpi,rupi ^S{cvpi)WNifjxapp^ > MS{cypi)p 

V Mu Hug 

'y / ^ ^ ^ f^■ugj ,rupi^vpi,uqjkQPtlsjgj ,rupi' V UjP, Z, fc 

U—1q—1j—1 

( 20 ) 

Constraints (|6]l and (|7]l with “if dyp > 0” removed; 
Constraints (|9]l, ( EOl l. (El, (El, (El, and (El. 


As an example, we consider the WPAN of six nodes showed 
in Fig. Etb). For the flow from A to B with 18 packets, we 
select three paths, path 1, path 2, and path 3, which have 
been illustrated in Fig. Elb). In this example, we assume for 
any two nonadjacent links, {v,p,i) and {u,qj), fsugjMvpi 
is equal to 0. Then we solve the problem of (P2) using the 
branch-and-bound method. The traffic distribution scheme is 
to send 3 packets through path 1, 9 packets through path 2, 
and 6 packets through path 3. The transmission phase of the 
schedule has 9 time slots, and is already illustrated in Fig.Efa). 
Compared with the single hop transmission scheme, MPMH 
reduces the number of time slots for transmission by 50%. 

The formulated mixed integer linear program (MILP) prob¬ 
lem (P2) is NP-hard. The number of decision variables 
is 0{{VPmaxPlmaxYK), and the number of constraints is 
0{{VPmaxHmaxYK), where Pmax Is the maximum number 
of paths of flows, and Hmax is the maximum number of 
hops of paths. Using the branch-and-bound algorithm to solve 
this problem will take significantly long computation time, 
e.g., several minutes for a 5-node network in ED, which 
is unacceptable for practical mmWave WPANs where the 
duration of one time slot is only a few microseconds ED- 
Therefore, to implement MPMH scheduling in a practical 
mmWave WPAN, heuristic algorithms with low computational 
complexity are needed. 


V. The MPMH Scheme 

To solve problem (P2), which boosts the spatial reuse via 
multiple paths multi-hop transmissions, firstly we should select 
suitable transmission paths. Then since the traffic distributed 
on these paths has a big impact on the efficiency of spa¬ 
tial reuse, traffic of each flow should be distributed on its 
transmission paths efficiently. Finally, transmission scheduling 
is needed to accommodate the traffic of all flows with a 
minimum number of time slots by fully exploiting spatial 
reuse. Following the above ideas, in this section, we propose 
a multiple paths multi-hop transmission scheme (MPMH), 
which consists of three heuristic algorithms respectively for the 
path selection, traffic distribution, and transmission scheduling. 
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A. Path Selection 

We set the maximum number of hops for each selected 
path, and denote it by Hmax- The path selection algorithm 
first finds out all the possible paths from the sender to 
the receiver with the number of hops less than or equal 
to Hmax, denoted by Pc{v). To exploit higher transmission 
ability of each path, each hop of these paths should have the 
same or higher channel quality than the direct path of flow 
V. Besides, each path should have no loop. Besides, if Cy is 
equal to 0, i.e., the direct path of flow v is blocked, flow v will 
not be transmitted through its direct path. Then the algorithm 
examines each path in Pc{v) in order of non-increasing lowest 
transmission rate on each path. The set of selected paths from 
Pc{v) is denoted by Psiv). Since the lowest transmission rate 
on each path determines its transmission ability, we should 
first select the paths with high transmission ability into Ps{v). 
Besides, the paths in Ps{v) should have no common hop to 
avoid degradation of efficiency. To maximize spatial reuse, 
the hops with the lowest transmission rate in Ps{v) should 
be nonadjacent to enable concurrent transmissions. Therefore, 
|Ps(i’)| is bounded by [n/2j with a network of n nodes. 

We denote the set of flows that will be transmitted through 
multiple paths as Fmpmh, which is selected according to (IHi. 
Then, the path selection algorithm should select suitable paths 
for each flow v in Fmpmh- We denote the sender of flow v as 
Sv, and the receiver of flow v as Vy. We also denote the first 
node and the last node of path p as fp and Ip respectively. For 
link from Ip to i, we denote its transmission rate by ci^i. The 
set of possible paths from Sy is denoted as P{v), and P{v) 
is initialized to the vertex of Sy. For each path p, we denote 
the lowest transmission rate on it as ci{p), and the hop with 
the lowest transmission rate as hi{p). We denote the set of 
lowest transmission rate hops of paths in Ps{v) is denoted as 
Hi{Ps{v)). 

The pseudo-code of the path selection algorithm is presented 
in Algorithm[T] Lines 1-18 And out all the possible paths with 
the number of hops less than or equal to Hmax, denoted by 
Pc{v). For each path in P{v), the algorithm extends this path 
to generate new paths with no loop, as indicated by lines 3-11. 
In line 5, the new hop extended from p should have a higher 
or same transmission rate than that of the direct path of flow v. 
After generating new paths from P{v), the old paths in P{v) 
are removed as in line 10, and P{v) is updated by the new 
set of paths P-new, as in line 12. For each path p in P{v), if 
its last node is the receiver r„, then this path will be added 
to Pc{v), and removed from P{v), as indicated by lines 14- 
16. Lines 19-28 select the eventual set of transmission paths, 
Ps{v), from Pc{v). The algorithm first obtains the path with 
the largest transmission ability, i.e., the lowest transmission 
rate on each path, as indicated by line 20. Then if this path 
has no common hop with the paths already in Ps{v), and the 
hop on this path with the lowest transmission rate, hi {p), is not 
adjacent to the hops in Hi{Ps{v)), this path will be selected 
into Ps{v), as indicated by lines 21-25. This path is removed 
from Pc{v) in line 26. When there is no path in Pc{v), the 
algorithm outputs Ps{v). 

For the WPAN in Fig. [Ub), with Hmax set to 3, the path 


selection algorithm obtains three paths for flow A ^ B, i.e., 
path 1, path 2, and path 3 already illustrated in Fig. [Hb). 


Algorithm 1 The Path Selection Algorithm 

Initialization: 

P{v) = {s„}; Pc{v) = 0; Psiv) = 0; iT;(Ps(u)) = 0; 
h=0; 

Iteration: 

1 : while (|P(t^)| > 0 and h < Hmax) do 

2: h=h-t-l, P-new — 05 

3: for each p € P(u) do 

4: for each node i with link Ip ^ i unblocked do 

5: if (Qj ^2 ^ Cy and i IS not on p) then 

6: Generate a new path p* by extending p to i; 

Pnew — Pnew G p , 

8: end if 

9: end for 

10: Piv) = Piv) — p; 

11 : end for 

12: Piv) — Pnew', 

13: for each p € Piv) do 

14: if {Ip == Vy) then 

15: Pciv) = Pciv) Up; Piv) = Piv) — p\ 

16 : end if 

17: end for 

18 : end while 

19: while (|Pc('f)| > 0) do 

20 : Obtain the path with the largest ciip), p G Pdv); 

21 : if (p have no common hop with paths in Psiv)) then 

22 : if (hiip) is not adjacent to the hops in HiiPsiv))) 

then 

23: Psiv) = Psiv) U p; HliPsiv)) = HliPsiv)) U 

hlip); 

24: end if 

25: end if 

26: Pciv) = Pciv) - p; 

21: end while 

28: Output Psiv). 


B. Traffic Distribution 

With the results of path selection, we propose a traffic 
distribution algorithm to distribute the traffic of a flow on 
the selected multiple paths, i.e., the traffic demand of flow 
V should be distributed on the set of selected paths, Psiv). 
Since the link with the lowest transmission rate on a path 
determines the transmission ability of the path, we should 
let the links with lowest transmission rates on different paths 
transmit concurrently as much as possible, and we also need to 
maximize the utilization of time slots in a pairing by making 
as many links as possible to transmit concurrently in each time 
slot of this pairing. Therefore, the traffic distribution algorithm 
is the proportional distribution of the traffic according to the 
lowest transmission rate on each path. For example, for the 
three paths in Fig. [lib), since the lowest transmission rates of 
path 1, path 2, and path 3 are 1, 3, and 2 respectively, the 








traffic from A to B is distributed on path 1, path 2, and path 
3 according to the proportion of 1:3:2. 

C. Transmission Scheduling 

After path selection and traffic distribution, the transmission 
scheduling algorithm computes a near-optimal schedule to 
accommodate the traffic demand of flows. Since any two 
adjacent links cannot be scheduled concurrently, it can also be 
inferred that the maximum number of links that can transmit 
in the same pairing is [n/2j ifTTI . Links that are scheduled 
in the same pairing can be represented by a directive graph, 
which is a matching mil. If we assign one of K colors to each 
pairing, and all the edges in the same pairing have identical 
color, each hop will have only one of the K colors since each 
hop is only scheduled in one pairing of the frame. Therefore, 
this process can also be modeled as an edge coloring problem 
mil, and the only difference is that the hops on the same path 
should be scheduled one by one from the first hop to the last 
hop. We denote the set of directional links of the pairing 
as iJ‘, and the set of vertices of the links in is denoted by 
L*. Thus, the problem of optimal scheduling is to find out the 
directive graph of each pairing to accommodate the traffic of 
all flows with a minimum of time slots. Our algorithm starts 
scheduling from the paths with the largest number of hops, 
and to maximize spatial reuse gain, the algorithm schedules 
as many links into each pairing as possible with the condition 
of concurrent transmission satisfied. Since the inherent order 
of hops in the same path, at most one hop can be scheduled 
in the same pairing, and the algorithm needs to schedule the 
preceding hops first. 

For each flow v, the set of its transmission paths is de¬ 
noted as Ps{v), which can be obtained by the path selection 
algorithm. For flows not in Fmpmh, its Ps{v) only has the 
direct path. From the traffic distribution algorithm, we obtains 
the traffic demand of each path in Ps{v). In the transmission 
scheduling, we denote the number of nodes as n, and the set 
of selected paths of all flows as Ps, which includes the paths 
in Ps{v) for each flow v. For each path p S Pg, we denote its 
number of hops as h{p). For each hop of each path, according 
to the traffic distributed on this path, we define the number of 
time slots to accommodate its traffic as the weight of this hop. 
For the hop link of path p € Pg, its weight is denoted as 
Wpi- We also denote the hop link of path p € Pg as {p, i). 
We denote the sender of link {p, i) by Spi, and the receiver by 
rpi. The set of hops in Pg is denoted by H. The hop number 
of the first hop that is not scheduled on path p is denoted as 
Fu{p)- In the pairing, the set of paths that are not visited 
yet is denoted by P*. 

The pseudo-code of the transmission scheduling algorithm 
is presented in Algorithm |2] First, we obtain the set of paths 
after path selection. From Pg, we can obtain the set of hops in 
Pg, FI. Since we should start scheduling from the first hop of 
each path, Fu (p) for each path p is set to 1. Then we iteratively 
schedule each hop of F[ into each pairing until all the hops in 
H are scheduled, as indicated in line 1. In each pairing, we 
first And out the unvisited paths with the maximum number 
of unscheduled hops as in line 6. Then among these paths. 


the first unscheduled hop with the minimum distance between 
its weight and current duration of this pairing is selected 
as the candidate hop of this pairing, as indicated in line 7. 
This step makes the numbers of time slots of hops in this 
pairing as close as possible, which can improve the utilization 
of time slots in this pairing as much as possible. Then the 
algorithm checks whether this candidate hop is adjacent to 
the hops already in this pairing since adjacent links cannot 
be scheduled concurrently, as indicated by line 8. If it is not, 
Arst the candidate hop will be added to this pairing, and then 
the concurrent transmission conditions of this pairing will be 
checked, as indicated by lines 9-16. If the conditions cannot 
be met, this candidate hop will be removed from this pairing, 
as indicated by lines 19-20. Otherwise, the duration of this 
pairing will updated to accommodate the traffic of this hop as 
in line 17. Since at most one hop of one path can be scheduled 
in one pairing, the path where this hop is will be removed from 
the set of the unvisit paths as in line 22. When the number of 
links in each pairing reaches [n/2j or there is no path in the 
unvisited path set, the algorithm will start scheduling for the 
next pairing as in line 5, and the scheduling for this pairing 
will be outputted as in line 24. 

Applying this algorithm to the example in section IIII-BI we 
obtain the schedule as follows: in the first pairing, link A —>■ D 
of path 3 transmits for one time slots; in the second pairing, 
link A —> 67 of path 2 and D ^ F of path 3 transmit for 
three time slots; in the third pairing, link A —> i? of path 1 
and C E of path 2 transmit for three time slots; in the 
fourth pairing, link F ^ B of path 3 transmits for one time 
slot; in the fifth pairing, link E ^ B transmits for two time 
slots. Thus, it needs 10 time slots to clear the traffic of A ^ B. 
As discussed before, the optimal solution needs 9 time slots. 
Thus our heuristic scheme obtains nearly the same scheduling 
solution. 

D. Computational Complexity and Control Overhead 

The path selection algorithm has the computational com¬ 
plexity of where n is the number of nodes in the 

network. For the traffic distribution algorithm, its computa¬ 
tional complexity is negligible. For the transmission schedul¬ 
ing algorithm, it has the complexity of 0(\Ps\rP), where Pg 
is the set of selected paths. Thus, the overall complexity of our 
scheme is -f \Ps\rP), which is a pseudo-polynomial 

time solution and suitable for the implementation in practical 
mmWave WPANs. 

With multi-path multi-hop transmissions, there will be in¬ 
creased control overhead since more scheduling information 
needs to be pushed to the nodes in the WPAN. However, the 
increase of the control overhead is small with respect to the 
Gbps transmission rates of mmWave links, and the influence 
on the performance is minimal. 

VI. Performance Evaluation 

In this section, we evaluate the performance of our proposed 
MPMH under various traffic patterns, and compare its perfor¬ 
mance with the optimal solution scheme and other existing 
protocols. 
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Algorithm 2 The Transmission Scheduling Algorithm 

Initialization: 

Input; the set of selected paths of all flows, 

Obtain the set of hops in Ps, denoted by if; 

Obtain the number of hops for each path p € Ps, h{p)-, 
Obtain the weight of each hop {p,i) G H, yUpP, 

Set Fu[p) = 1 for each p G Ps', f=0; 

Iteration: 

1 : while (|iT| > 0) do 
2: f=f+l; 

3: Set V* = 0, = 0, and = 0; 

4: Set P* with P^ = Ps; 

5: while (|P^| > 0 and |P*| < \n/2\) do 

6 : Get the set of unvisited paths with the largest number 

of unscheduled hops, Pmh', 

7: Get the hop (p, P„(p)) of path p G Pmh with the 

minimum abs((5* — Wpp^(^p-^)', 

8 : if (SpF„(p) ^ y* and ^ V*) then 

9 : P* =P‘u{(p,P,(p))}; 

10- y y ^ {'^pFu(p) ppFu (p)}’ 

11 : for each link {p,i) in H* do 

12 : Calculate the SINR of link {p,i), SINRpi 

13: if (SINRpi < MS{cpi)) then 

14: Go to line 19 

15: end if 

16 : end for 

17: (5* = max(b‘, tUpF„(p)), H = H - {p, P«(p)); 

18: Pu(p) = Fu{p) + 1; Go to line 21 

19: P‘=P‘-{(p,P„(p))}; 

20 ; y y i-^pFu{p)NpFuip)}•> 

21 : end if 

22: P‘ = P‘-p; 

23: end while 

24: Output H* and b*; 

25: end while 


A. Simulation Setup 

In the simulation, we consider a typical mm Wave WPAN of 
10 nodes. We assume that all the WNs and the PNC are uni¬ 
formly distributed in a square area with 8m x 8m. According 
to the distances between WNs, we set four transmission rates, 
2 Gbps, 4 Gbps, 6 Gbps, and 8 Gbps. There are 10 flows in the 
network. With e of the criterion in Section ITV-AI set to 0.0625, 
there is only one flow transmitted through multiple paths, and 
other flows are transmitted through the direct paths. The paths 
are selected according to the algorithm in Section [V-AI We set 
the size of data packets to 1000 bytes. We adopt the simulation 
parameters in Table II of ||20l, which is also summarized in 
Table U The duration of one time slot is set to 5 ps. With 
the transmission rate of 2 Gbps, a packet can be transmitted 
in one time slot. As in M, for the simulated network, the 
PNC can complete the polling of traffic or schedule pushing 
in one time slot. Generally, it takes only a few time slots for 
the PNC to complete path selection and schedule computation. 
To adapt to dynamical network states, the number of time slots 
of a frame is bounded by 1000. The simulation length is set to 


TABLE I 

Simulation Parameters 


Parameter 

Symbol 

Value 

PHY data rate 

R 

2 Gbps, 4 Gbps, 6 Gbps, 8 Gbps 

Propagation delay 

5p 

50ns 

Slot Duration 

Tslot 

5 ps 

PHY overhead 

TpHY 

250ns 

Short MAC frame Tx time 

TshFr 


Packet transmission time 

'Vpacket 

1000*8/R 

SIFS interval 

TsifS 

100ns 

ACK Tx time 

Tack 

TshFr 


5 X 10^ time slots. We also set the delay threshold to 2.5 x 10^, 
and will discard packets with delay larger than the threshold. 
Initially, each flow has a few packets generated randomly. The 
maximum number of hops for each selected path, Hmax, is 
set to 3. In the simulation, we assume nonadjacent links are 
able to be scheduled for concurrent transmissions. 

In the simulation, we set the following two kinds of traffic 
modes: 

1) Poisson Process: packets arrive at each flow following 
a poisson process with arrival rate A. The traffic load in this 
mode, denoted by Ti, is defined as 


Ti 


X X L xV 
R 


( 21 ) 


where L is the size of data packets, V is the number of flows, 
and R is set to 2 Gbps. 

2) Interrupted Poisson Process (IPP): packets arrive at 
each flow following an interrupted poisson process (IPP). The 
parameters of the interrupted poisson process are Ai, A 2 , pi 
and p 2 , and the arrival intervals of an IPP obey the second- 
order hyper-exponential distribution with a mean of 


E{X) 


El + El 

Ai A 2 


( 22 ) 


The interrupted Poisson process can be represented by an 
ON-OFF process. The ON duration and the OFF duration obey 
the negative exponential distribution of ri and ^ 2 , respectively. 
In the ON duration, packets arrive at each node following the 
Poisson process of arrival rate Xon_off, while no packet arrives 
in the OFF duration. Xon_off^ ^i, and r 2 can be inferred from 
Ai, A 2 , Pi, and p 2 as follows. 


^on_off — PlXl + P2X2, 


n 


PlP2{Xl — A2)^ 


A, 


on_off 


7-2 


Ai A2 

^on_off 


(23) 

(24) 

(25) 


Therefore, IPP traffic is typical bursty traffic. We define the 
traffic load Ti in this mode as; 


Ti = 


LxV 
E{X) X R' 


(26) 


We evaluate the system by the following four performance 
metrics; 

1) Average Transmission Delay: The average transmission 
delay of received packets from all flows, which is in units of 
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time slots. The average transmission delay does not include 
the time caused by steering antenna. 

2) Network Throughput: The number of successful trans¬ 
missions of all flows until the end of simulation. For a 
packet of one flow, if its delay is less than or equal to the 
threshold, it will be counted as a successful transmission. With 
constant simulation length and fixed packet size, total number 
of successful transmissions is a good indication to show the 
throughput performance. 

3) Average Flow Delay: The average transmission delay of 
the flow transmitted through multiple paths. 

4) Flow Throughput: The number of successful transmis¬ 
sions achieved by the flow transmitted through multiple paths. 

In order to show the advantages of MPMH, we compare 
MPMH with the following two protocols: 

1) FDMAC: The frame-based scheduling directional MAC 
protocol, and the core of FDMAC is the greedy coloring (GC) 
algorithm, which can compute near-optimal schedules with 
respect to the total transmission time with low complexity im. 
FDMAC is also a frame based protocol, and in GC, the traffic 
demand of flows is accommodated by iteratively scheduling 
each flow into each concurrent transmission pairing in non¬ 
increasing order according to traffic demand. To the best of 
our knowledge, FDMAC achieves the best performance among 
the existing protocols, such as MRDMAC ll20l and MDMAC 
(El, by fully exploiting spatial reuse in mmWave WPANs. 

2) FDMAC-UR: FDMAC with links’ differences in the 
transmission rate not considered, and the transmission rates 
of all links are set to 1 Gbps. 

B. Comparison with the Optimal Solution 

We first compare MPMH with the optimal solution of prob¬ 
lem (P2) by the branch-and-bound method. Since obtaining the 
optimal solutions takes extremely long time, we assume there 
is only one flow with the traffic to transmit, and the traffic of 
this flow needs to be transmitted through multiple paths. The 
traffic load is defined as in (l2Tl i and (|2^ with V equal to 1. 
In this case, the delay threshold is set to 3 x 10"^. 

The average flow delay and the flow throughput of MPMH 
and the optimal solution are plotted in Fig. |2] From the 
results, we can observe that the gap between MPMH and 
optimal solution is negligible under light load. For the average 
flow delay, the gap at the traffic load of 5 is only 7.9% of 
the average flow delay of MPMH. For the flow throughput, 
the gap is only 5.3% of the flow throughput of MPMH. 
Therefore, we have demonstrated that MPMH achieves near- 
optimal performance with low complexity. 



(a) Average flow delay (b) Flow throughput 

Fig. 2. Comparison between MPMH and the optimal solution under Poisson 
traffic. 


The average execution time of schedule computation for 
MPMH and the optimal solution is plotted in Fig. [3 From the 
results, we can observe that MPMH has much less execution 
time than the optimal solution, and thus it has much lower 
computational complexity. 



Fig. 3. The execution time of MPMH and the optimal solution under Poisson 
traffic. 

C. Comparison with Existing Protocols 

1) Delay: We then evaluate the average network delay of 
the three MAC protocols. In Fig. (H we plot it for different 
traffic loads. From the results, we can observe that with 
the increase of traffic load, the delay of the three protocols 
increases, and the delay of FDMAC-UR increases rapidly 
under light load. MPMH and FDMAC have similar delay 
performance under light load. Under light load, the arrived 
packets can be transmitted in a short time, and the frame 
length is short. Thus, the delay in this case is small. Under 
heavy load, MPMH outperforms FDMAC and FDMAC-UR. 
With the traffic load from 4 to 7, compared with FDMAC, 
MPMH decreases the average transmission delay by about 
75.74% under Poisson traffic and 86.54% under IPP traffic 
on average. For FDMAC-UR, since the actual transmission 
ability of links is not exploited fully, it has much higher delay 
compared with MPMH and FDMAC. These phenomena can 
be explained as follows. In FDMAC and FDMAC-UR, flows 
cannot be transmitted through multiple paths of multiple hops. 
Thus, flows with low channel quality on their direct paths 
will occupy a large number of time slots in the schedule, 
which increases the frame length significantly and also delay 
of packets. Since the frame length is bounded by 1000 time 
slots, the packets that cannot be transmitted are re-scheduled 
in the next frame. Thus, with the increase of the traffic load, 
the system enters saturation, and the curves become flat and 
show a concave form. 



Fig. 4. Average transmission delay of the three MAC protocols under 
different traffic loads. 
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As regard to the average flow delay, we present the results 
in Fig. 13 We can observe that the delay curves of MPMH and 
FDMAC begin to diverge at the traffic load of 3. Compared 
with FDMAC, MPMH decreases the delay by about 74.31% 
under Poisson traffic and 74.29% under IPP traffic on average 
with the traffic load from 4 to 7, respectively. 



(a) Poisson traffic (b) IPP traffic 

Fig. 5. Average flow delay of the thi'ee MAC protocols under different traffic 
loads. 


2) Throughput: The network throughput achieved by the 
three protocols is plotted in Fig. |6] We can observe that 
MPMH has the maximum throughput in all cases, and the 
gap between MPMH and FDMAC is more significant under 
heavy load. Under light load, the delay is small, and all the 
arrived packets can be transmitted successfully for MPMH 
and FDMAC. Thus, the throughput of MPMH and FDMAC 
increase linearly under light load. FDMAC stops increasing at 
the traffic load of 4, while MPMH stops increasing at the traffic 
load of 6 since the network tends to saturation. For FDMAC- 
UR, its throughput is poor, and tends to be saturated at the 
traffic load of 2. Compared with FDMAC, MPMH increases 
the network throughput with the traffic load from 5 to 10 
on average by about 54.37% under Poisson traffic and about 
50.58% under IPP traffic, respectively. When traffic load is 10, 
MPMH outperforms FDMAC by about 80.2%. The reason for 
MPMH to outperform FDMAC is that in FDMAC, flows with 
low channel quality on their direct paths cannot be transmitted 
through multiple multi-hop paths as in MPMH, and occupy 
a large number of time slots in the schedule, which are 
underutilized for concurrent transmissions. Under light traffic 
load, the throughput of MPMH and FDMAC is determined 
by the traffic arriving at each node. Under heavy traffic load, 
however, the advantages of MPMH over FDMAC will stand 
out, and MPMH will outperform FDMAC significantly. On the 
other hand, with the increase of traffic load, delays of packets 
increase, and a considerable number of packets are discarded 
due to their delays exceeding the threshold, which leads to 
a decrease in FDMAC throughput from the traffic load of 7 
and the bigger gap between MPMH and FDMAC at the traffic 
load of 10. 

On the other hand, the flow throughput of the three protocols 
is presented in Fig. [T] We can observe that the tendencies of 
the curves are similar to those in Fig. |6] With the traffic load 
from 5 to 10, MPMH outperforms FDMAC by about 52.14% 
under Poisson traffic and 47.66% under IPP traffic on average, 
respectively. The result show that MPMH greatly increases the 
throughput of the flow of low channel quality compared to 
FDMAC, especially under heavy traffic loads. 




(a) Poisson traffic (b) IPP traffic 

Fig. 6. Network throughput of the three MAC protocols under different 
traffic loads. 



(a) Poisson traffic (b) IPP traffic 

Fig. 7. Flow throughput of the three MAC protocols under different traffic 
loads. 


D. Performance under different Hmax 

To investigate the impact of the choice of Hmax on network 
performance, we plot the average transmission delay and 
network throughput of MPMH with Hmax equal to 2, 3, and 
4, respectively, in Fig. 0 The traffic mode is the Poisson 
traffic. From the results, we can observe that MPMH with 
Hmax equal to 3 achieves the best performance in terms of 
delay and throughput. However, the gap between MPMH with 
different Hmax is small. For MPMH with Hmax equal to 2, 
the advantage of multi-path multi-hop transmissions is limited 
by the number of hops on paths. For MPMH with Hmax equal 
to 4, more hops will lead to larger delay, and with the delays 
of more packets exceed the threshold, the network throughput 
also degrades. Therefore, the maximum number of hops on 
each path, Hmax, should be optimized to achieve an optimal 
network performance in practice. 



(a) Average Transmission Delay (b) Network Throughput 

Fig. 8. Delay and throughput of MPMFl with different Hmax under Poisson 
traffic. 

We also plot the flow delay and throughput of MPMH with 
Hmax equal to 2, 3, and 4, respectively, in Fig. |9l We can 
observe that the results are similar to those in Fig. [3 and 
MPMH with Hmax equal to 3 has the best performance in 
terms of flow delay and throughput. With a small Hmax as 
2 , the number of paths for each flow is limited, and the 
concurrent transmissions among paths cannot be exploited 
fully to improve flow delay and throughput performance. With 
a large Hmax as 4, transmissions through paths with more hops 
lead to larger delay, and there are more packets discarded since 
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their delays exceed the threshold. In practice, Hmax should be 
selected according to the actual network conditions to optimize 
network performance. 




(a) Average Flow Delay (b) Flow Throughput 

Fig. 9. Flow Delay and throughput of MPMH with different Hmax under 
Poisson traffic. 

VII. Conclusion 

In this paper, we proposed MPMH for mmWave WPANs in 
the 60 GHz band, which boosts the potential of spatial reuse 
by transmitting the traffic of flows with low channel quality on 
their direct paths or with high traffic demand through multiple 
multi-hop paths to improve throughput and to reduce latency 
for both these flows and the network. Extensive simulations 
demonstrate that MPMH achieves near-optimal performance 
compared with the optimal solution. Compared with FDMAC, 
MPMH increases flow and network throughput by about 
50% and 52.48% on average, respectively. Performance under 
different maximum number of hops indicates the maximum 
number of hops should be selected according to the actual 
network conditions to optimize network performance. 

In the future work, we will investigate the accurate and 
reasonable modeling of NLOS links and try to incorporate 
NLOS transmission into our scheme. In MPMH, the choice of 
£ controls the number of flows transmitted through multiple 
paths, and we will further optimize e according to the actual 
network conditions such as the channel transmission rate 
distribution of links and the traffic demand distribution of 
flows. Although the fairness among flows is improved by 
the multi-path multi-hop transmission, we will investigate the 
mechanisms to improve the fairness of our scheme further. 
Besides, we will also evaluate the energy consumption and 
energy efficiency of MPMH. 

Appendix A 

Notation in Problem Formulation 

In order to facilitate the reader to understand the notations 
in Problem Formulation, we list the notations in Table |II] as 
follows. 


Appendix B 
Notation in MPMH 

To facilitate the understanding of MPMH, we also list the 
notations in Table |III] as follows. 


TABLE II 

Notation in Problem Formulation 


Symbol 

Description 

dy 

The traffic demand of flow v 


The number of paths of flow v 

Hyp 

The number of hops of the path of flow v 

dyp 

The traffic demand distributed to the path 

{v,p,i) 

The hop link of the p*^ path of flow v 

t^vpi 

The transmission rate of (i»,p, 2 ) 

dvpi,uqj 

A binary variable to indicate 
if and {u,q,j) are adjacent 

K 

Number of pairings in a schedule 

5>^ 

Number of time slots of the pairing 


A binary variable to indicate whether link {v,p,i) 

vpi 

is scheduled in the pairing 

^vpi 

Sender of link (i;, p, i) 

't' vpi 

Receiver of link (v,p,i) 

f ^uqj ) ^vpi 

A binary variable to indicate whether Suqj and r^pi 
direct their beams towards each other 


The distance between Sypi and r^pi 



TABLE III 

Notation in MPMH 

Symbol 

Description 

Hmax 

The maximum number of hops for each selected path 

Sy 

The sender of flow v 

Tv 

The receiver of flow v 

P{v) 

Set of all possible paths from Sy 

Hnew 

Set of generated new paths from P{v) 

Pciv) 

Set of all possible paths from Sy to Py 

Hmpmh 

Set of flows transmitted through multiple paths 

fp 

The first node of path p 

Ip 

The last node of path p 

Clip) 

The lowest transmission rate on path p 

hlip) 

The hop on path p with the lowest transmission rate 

Psiv) 

The set of selected paths 

HliPsiv)) 

The set of lowest transmission rate hops of paths in Psiy) 

S* 

Number of time slots of the pairing 


Set of directional links in the pairing 

u* 

Set of vertices of the links in 

n 

The number of nodes 

Ps 

Set of selected paths of all flows 

h{p) 

Number of hops on p 

iP,i) 

The hop link of path p 

^pi 

The sender of link (p, i) 

t'pi 

The receiver of link (p, i) 

'^pi 

The weight of link (p, i) 

H 

The set of hops in Ps 

Fu[p) 

The hop number of the first unscheduled hop on path p 

pt 

The set of unvisited paths 

Hmh 

Set of unvisited paths with the lai'gest 
number of unscheduled hops 
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